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ABSTRACT 

The low order moments, such as the bulk flow and shear, of the large scale peculiar 
velocity field are sensitive probes of the matter density fluctuations on very large 
scales. In practice, however, peculiar velocity surveys are usually sparse and noisy, 
which can lead to the aliasing of small scale power into what is meant to be a probe 
of the largest scales. Previously, we developed an optimal "minimum variance" (MV) 
weighting scheme, designed to overcome this problem by minimizing the difference 
between the measured bulk flow (BF) and that which would be measured by an ideal 
survey. Here we extend this MV analysis to include the shear and octupole moments, 
which are designed to have almost no correlations between them so that they are 
virtually orthogonal. We apply this MV analysis to a compilation of all major peculiar 
velocity surveys, consisting of 4536 measurements. Our estimate of the BF on scales of 
~ lOO/i^Mpc has a magnitude of \v\ = 416 ±78 km/s towards Galactic / = 282° ±11° 
and b = 6° ±6°. This result is in disagreement with ACDM with WMAP5 cosmological 
parameters at a high confidence level, but is in good agreement with our previous 
MV result without an orthogonality constraint, showing that the shear and octupole 
moments did not contaminate the previous BF measurement. The shear and octupole 
moments are consistent with WMAP5 power spectrum, although the measurement 
noise is larger for these moments than for the BF. The relatively low shear moments 
suggest that the sources responsible for the BF are at large distances. 



Subject headings: cosmology: distance scales - cosmology: 
large scale structure of the universe - cosmology: observation 
- cosmology: theory - galaxies: kinematics and dynamics - 
galaxies: statistics 



1 INTRODUCTION 

Large-scale structure formation is assumed to arise from 
small Gaussian ini tial fluctuations amplified by gravi- 
tatio nal instability (B ardeen et alj 1 19861 : lEisenstein fe Hul 
1998). This basic framework is strongly supported by the 
consistency between the Cosm ic Microwave Backg round 
(CMB) angular power spectra (|Dunklev et alj 120091 here- 
after WMAP5) observed at high-redshift, and l arge-scale 
structure (LSS) data from gravitational lensing (|Fu et alj 



20081) and galax y power spectra (lEisenstein et al.1 120051; 
Cole et al. 20051) and its bispectrum (IScoccimarro et al.l 



200 ll : iFeldman et al.ll200ll : IVerde et alj|2002h , measured at 



low redshift. 

However, on the largest scales, the comparisons are usu- 
ally between the matter density fluctuations measured by 
the CMB and the fluctuations in galaxy number density, 
and so are susceptible to uncertainties in the relationship 
between mass and light ( "biasing" ) . While dark matter can 
be observed directly at relatively low redshift via gravita- 
tional lensing, at pre sent this techn ique only just reaches 
into the linear regime (|Fu et al.ll200S ) . The Integrated Sachs- 
Wolfe (ISW) effect (|Sachs fc Wolfe! 1 19671 ) probes the evo- 
lution of the dark matter potential on large scales via a 
cross-correlaton between galaxies and the CMB. Recent 
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ISW compilatio ns suggest a stronger signal than expected 
(|Ho et alj 12008V Perhaps the most promising method for 
probing the low-redshift DM power spectrum on very large 
scales (< 100/i _1 Mpc, where h is the Hubble constant in 
units of 100 km s" 1 Mpc - 1 ) is via the peculiar v elocity 
field of galaxies and clusters ((Strauss fc Willickll 19951 ). 

Peculiar velocities can be measured statistically 
through redshift-space distortions of gala xy power spectra 
jPeacock et al.ll200ll ; iTegmark et~ai]|2004r i , or more directly 
by measuring distances to individual galaxies or clusters. 
Recent peculiar velocit y studies, using vari o us techniques 
dPike fc Hudson! 120051; iPark fe Park! 120061; ISarkar et al l 



20071; IWatkins fe Feldman! 120071 ; iFeldman fe Watkins! 120081 ; 
Watkins et "all I2009T) . all suggest that different peculiar 
velocity surveys, using different distance estimators, are 
all consistent with sampling the same underlying peculiar 
velocity field. Moreover, recent peculiar v elocity surveys 
dMasters et al.l 120061 ; ISpringob etafl 120071 ; ISpringob etakl 
2009) are deeper, denser, and more reliable than ever be- 
fore. As surveys have gotten larger, our understanding of 
the distance indicators needed to extract the peculiar veloc- 
ities, and to control their systematic errors, has also im- 
proved. Finally, new analytic techniques h ave allowed us 
to be t ter extract information from s urveys jFeldman et alj 
20031; iRadburn-Smith et al l 12004 IPike fe Hudson! 120051; 



Sarkar et ai1l2007l ; I Watkins fc Feldman! 120071 ; lLavaux et"al] 
2010). 



On small scales (<J 20 h 1 Mpc), peculiar velocities yield 
results consistent with the WMAP5+B AO+SN ACDM 
mode l value: (n m /0.3)°- 55 a 8 = 0.77 ±0.035 jKomatsu et al l 



2009). For example, on these scales, IPike fc Hudson! |2005h 
found (n m /0.3)°- 55 o- 8 = 0.80 ± 0.05 from a compar- 
ison of density (galaxy redshift) and peculiar veloc- 
ity surveys. A statis tical analysis of pairwise velocities 
(jFeldman et al. 20031) yields, afte r correction for non- 
linearities ( Juszkiewicz et ail 120101 1 , a slightly higher re- 
sult for as: 1-02 ± 0.18 , but which is sti l l con sistent 
with lKomatsu et all |2009l ). lAbate fc Erdogdul (I2009T ) found 
(fi m /0.25) a55 o-g = 0.901^16 from a correlation function 
analysis of SFI++ (Spiral Field I band) over a range of scales 
^ 25 h^Mpc. 

On the very largest scales, however, the agreement 
is less clear. Previous work has suggested that much of 
the peculiar velocity of the Local Group (LG) and nearby 
galaxies is generated by gravitational sources on very large 
scales, at a level which may be in excess of expectations 
from the standard ACDM model. This evidence comes di- 
rectly from estimates of th e bulk flow (BF) on la rge scales 
(IFeldman fc Watkins! 120081 ; IWatkins et all I2009T I. It also 
comes indirectly from galaxy redshift surveys which allow 
one to estimate how much of the LG's motion is generated 
locally: the measurement of a low amplitude, misaligned 
gravity dipole from to masses within J> 60 — 100 /i _1 Mpc 
then requires that a large component of the LG's motion is 
generated external ly, on larger scales (|Pike fc Hudsonl feoOS; 
lLavaux et aT1l2010l ')- 

Many previous studies of peculiar velocity data have 
been based on maximum-likelihood analysis of the entir e 
velocity field |Zaroubi et al.ll200ll ; lAbate fc Erdogdul 120091 ). 



While such studies have the advantage that they use all of 
the information present, there are disadvantages that arise 
because all scales are analyzed simultaneously. In particu- 
lar they are sensitive to treatment of quasi- and non-linear 
regimes, and to the details of the assumed peculiar veloc- 
ity errors. Moreover, such studies typically assume a given 
power spectrum shape a priori. In this work, we isolate the 
largest scales. 

In a recent paper (Watkins et al. 2009, hereafter Paper 
I), we developed the "minimum variance" (MV) moments 
that were designed to estimate the BF on a particular scale 
with minimal sensitivity to small-scale power. Paper I also 
showed that the BF's from independent peculiar velocity 
surveys were consistent with each other. In this paper, we 
extend the formalism to include the next higher elements 
in the expansion, namely the shear and octupole moments. 
These higher-order moments contain information about the 
power spectrum on scales that are large, but not as large as 
that probed by the BF. The primary goal of this paper is to 
assess whether these moments, like the BF, have higher am- 
plitude than expected than expected in ACDM. Moreover, 
they also all ow us to e x tract cosmographical information: 
for example, iLilie et all (| 19861 ) first used the existence of a 
BF and a shear in the very nearby Universe to determine 
the approximate distance to the Great Attractor. 

It is important to note that in our analysis, since we are 
fitting to moments of an idealized survey, our model does not 
change with additional moments, that is, we are estimating 
individual moments rather than fit a model. In principle, 
then, our estimates of bulk flow from Paper I should not 
change in the analysis presented here. However, in practice, 
our use of an orthogonality constraint as described in Sec- 
tion [2] below will lead to small changes in our bulk flow 
estimates as we include higher moments. 

In Section[2]we introduce the MV weights and the nine- 
teen MV moments for the BF, shear and octupole. In Sec- 
tion [3] we discuss the peculiar velocity catalogues analyzed 
here. In Section [4] we present the moment amplitudes and 
compare these with expectations from cosmological models. 
We discuss our results in Section[5]and conclude in Section[6] 



2 MOMENTS OF THE VELOCITY FIELD 

While large-scale flows are still in the linear regime, the 
statistics of individual galaxy or cluster peculiar velocities 
Sn are not well described by linear theory due to the exis- 
tence of nonlinear flows on small scales. This problem may 
be solved by decomposing the velocity field into components, 
each of which reflects motions on a particular range of scales. 
Since the statistics of moments associated with large-scale 
motions can be treated using linear theory, the amplitudes 
of these moments can be used to put direct constraints on 
cosmological parameters. 

The most commonly used decompositi on of the ve- 
locity field is a T aylor series expansion (|Kaiserl 1 19881 ; 
Ijaffe fc Kaiserl fl995l ). where the components of the veloc- 
ity field are written as 



Vi(r) = Ui + Uijrj + Uij k rjr k + . 



(1) 
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The three zeroth order constants, Ui, are the components of 
the oft-discussed "bulk flow" (BF). The first order "shear" 
tensor, U%j, is symmetric if the velocity field is assumed to 
be curl-free, as it must be if motions are caused by gravi- 
tational instability. Thus the condition that L7y = Uji gives 
six independent shear components. The second order ten- 
sor, Uijk, which we will call the "octupole" tensor, must 
also be symmetric for the same reason, so that there are 
ten independent octupole components, giving a total of 19 
independent components for a second-order expansion. 

This expansion is used to decompose the velocity field 
in a particular volume, typically that occupied by the galax- 
ies in a peculiar velocity survey. Under these circumstances, 
the BF moments probe scales larger than the diameter of 
the volume, with each subsequent order probing smaller and 
smaller scales. However, peculiar velocity surveys typically 
have complicated geometries, so that mom ent amplitudes 
are not comparable be tween surveys (see IWatkins et alJ 
120091 : ISarkar et alj|2007h . Indeed, even the interpretation of 
these moments can be difficult due not only to the compli- 
cated distribution of survey objects but also to the varying 
measurement errors associated with each object. 

Interpretation of peculiar velocity data would be much 
more straightforward if the data were close to an "ideal sur- 
vey": an infinitely large spherically symmetric survey with 
no measurement errors and with a gaussian radial distri- 
bution function f(r) oc e~ T ' 2R i , where the parameter Ri 
designates the depth of the survey. The velocity moments ob- 
tained in this way correspond to a known scale, are straight- 
forward to interpret, and are comparable between surveys. 
While we do not have such ideal surveys, we have some flex- 
ibility to force our actual surveys to match such an ideal 
survey as closely as possible. This flexibility is in the form 
of a weight for each peculiar velocity datum, which we are 
free to adjust in an optimal way. 

Note that this approach of adjusting the weights to 
match a given "ideal "geometry is very different from the 
standard maximum likelihood (MLE) weights, which min- 
imize the measurement error but do not account for the 
geometry of the survey. 

Before we discuss the estimation of the velocity mo- 
ments using velocity surveys, we must address a problem 
with the Taylor expansion that arises at second-order and 
beyond. While the zeroth and first order moments are or- 
thogonal, there is significant overlap between the BF and 
the octupole moments. In particular, a pure octupole flow 
in a given volume V, Vi = UijkriVk, contains a net BF in 
that volume given by J y Uij^TiTk d 3 r. This leads to a strong 
correlation between the BF and octupole moments. 

We can solve this problem by modifying the definition 
of the octupole moments. We rewrite our expansion of the 
velocity field in a volume V as 

Vi(r) = Ui + UijTj + Uijk {rjr k - A jk ) + ... (2) 
where the constants Ajk are given by 

Aj h = / TjTk d 3 r (3) 
Jv 

For a spherically symmetric volume, as we use here, only 



the diagonal elements An, A22, and A33 are nonzero. When 
the octupole moments are defined in this way, the 19 mo- 
ments to second order are orthogonal and have no overlap. 
A similar procedure can be carried out to remove overlap 
between the third order moments and the shear, however, 
current surveys are not sensitive to higher moments than the 
octupole and thus in this paper we will focus on the second 
order expansion only. 

Now that we have defined the velocity moments we are 
interested in, we turn to how these moments can be esti- 
mated using peculiar velocity data. First, it is only possible 
to measure the line-of-sight peculiar velocity s = v • f. Our 
expansion for the peculiar velocity field thus translates into 
an expansion for the line-of-sight velocity field s(r), which 
can be written as 

s(r) = Uih + Uijrfifj + Uijk {r 2 fifjf k - Ajkh) + • ■• (4) 

For simplicity, we follow [jaffe fc Kaiserl l|l995l ) and write the 
second-order expansion in the form of a 19 component vector 
of moment amplitudes, 
19 

S (r) = ^C/ pffp (r) (5) 

p=i 

where U p are the 19 moment amplitudes given by 

Up {Ux> Uy , Uz , UxX? Uyy , Uzz> U xy i Uyz? U ZX (6) 

U XXX , Uyyy , UzZZi U xxy i Uyyzj UzzX> U xyy ■> 

Uyzzi U zxx i Uxyz\ 

and the mode functions are given by 

Sj>( r ) = {fx,fy,fz,rfx,rfl,rf 2 z,2rf x r y ,2rfyfz, (7) 

37" TyTz ft^yyTzi 3t T z T x ^zz^xi 
37" TyTx ^yyfx-i 3/* T ' z Ty -A. Z z7*y, 

3r 2 f x fz - Axxf z ,6r 2 fxfyf z } 

where we have used the fact that only the diagonal elements 
of Aij are nonzero. 

We first consider an idealized survey, consisting of posi- 
tions r n and exact line-of-sight velocities s n for a spherically 
symmetric distribution of N objects with a given distribu- 
tion function f(r). In this case, the ideal velocity moments 
Up are just the projections of the velocities onto the mode 
functions, 

1 N ° 

Up = jj-^2g P (rn)s n (8) 

° n=l 

Thus the moment amplitudes take the form of linear com- 
binations of the velocities X/n w p." Sn w ith the numerical 
values of the MV weights given by 

™'p,n = 9 P (rn)/N (9) 

We note here that Eqn. ((8| applies only for moments that 
have no overlap, as is true for the moments defined in Eqn. 
0. The expression for moment amplitudes when overlap- 
ping moments are used is somewhat more complicated. 
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An actual peculiar velocity survey consists of N objects 
with positions r n and measured line-of-sight velocities S„ 
with uncertainties a n ■ The measured velocities are assumed 
to have the form S„ = s n + fin, where S n is drawn from a 
Gaussian distribution of variance a n + a 2 . Here <r* is the 
velocity noise, which accounts for small-scale motions not 
included in the measured moment. In the peculiar velocity 
literature, the approach often taken is one where one fits a 
flow model (whether parametric flow model, e.g. an expan- 
sion over many components, or using a "template" based on 
e.g. the gravity of a galaxy density field). In such cases, it 
is important that the flow model is complete on all scales 
and that the noise estimates and % 2 values are reasonable. 
That is not the approach taken here: instead we are esti- 
mating large-scale moments, which are expected to be close 
to orthogonal. Hence the value of a, affects our moments 
only very weakly: it modifies the measurement noise and 
hence the weights. While it has the strongest effect on nearby 
galaxies for which the velocity errors are smallest, these same 
galaxies are strongly downweighted by our Ri = 50 ft, - Mpc 
MV weighting scheme. Consequently, changing the value of 
<r* does not alter any of our moments significantly. 

Given an idealized survey with velocity moments Up, 
we wish to determine the weights w p , n such that the linear 
combinations 



n 



To find the weights that minimize this expression, we 
take the derivative with respect to w p , n , set it equal to zero: 

(S n S m ) +J2 X Pi9<i{rn) = (14) 

m q 

and solve for the weights, 

{SmU p )-^y ' \ pq g q (r m ) I (15) 



Up,n — ^ ^ G n m 



9 / 



where G is the covariance matrix of the individual measured 
velocities, G nm = (S n S m ). The values of the Lagrange mul- 
tipliers can be found by plugging Eqn. (|15|) into Eqn. (|11[) 
and solving for \ pq , 



x pi — 



A/^ 1 ^ Gnm (S m Ul)g q (r n 



where the matrix M is given by 
M pq = ^ ^ Gn m g p (r„)g q (r m ) 



(16) 



(17) 



n=l 



(10) 



give the best possible estimates of the U p . Following our 
previous work (Paper I), we calculate the weights by mini- 
mizing the average variance ((U p — u p ) 2 ). When considering 
only the BF moments in Paper I, we included additional 
constraints that ensured that the estimators would give the 
correct amplitude for a pure BF velocity field. Here we im- 
plement a more general set of constraints that is applicable 
to higher-order velocity moments. 

Suppose that the flow field consisted only of BF, shear, 
and octupole moments, so that the line-of-sight velocities at 
positions r„ took the form s„ = ^ p U p g p (r n ). In order for 
the estimators to give the correct amplitudes for the velocity 
moments on average for this flow field, i.e. (u p ) — U p , we 
require that 



ig q (r n ) = 5 P 



(11) 



This set of constraints can be implemented by using La- 
grange multipliers. Thus we seek to minimize the quantity 



{(Up - u p ) 2 ) + ^ \ pq I w p , n g q (v n ) - S pq ) 

q \ n / 



(12) 



Expanding out the first term, plugging in the expression for 
Up from Eqn. H10|) . we can write this expression in terms of 
the weights w p , n , 

{U 2 ) — 2w p , n {S n U p ) + w p ,„w p , m {S„S m ) (13) 



Eqn. (|15l) gives us a formula for calculating the MV 
weights in terms of the covariance matrix Gnm = {SnSm} 
and the correlation {S m U p ), both of which can be calcu- 
lated given a power spectrum model. Using the fact that 
S n = s n + 5 n as described above and that s n and S n are 
independent, we can write the covariance matrix as 

Gmn = (s m Sn) + fining + v'l) (18) 

= (f m -v(r m ) ?„ • v(r n )) + 8 m n{ai + &n)- 

In linear theory the first term can be expressed as an integral 
over the density power spectrum P(k), 



<f m ■ v(r 

m) I*7i " v(r„)> 



2^2 



dk P(k)fmn(k) (19) 



where the function / m „(fc) is the angle averaged window 
function 



/mn(fc) = J (?m ■ t) (? n ■ k) 

x exp (ik k ■ (r m - r n )j 



(20) 



The correlation {S m U p ) is calculated in a similar way. 
We generate an ideal survey by selecting N random posi- 
tions r' n , with the desired radial distribution function. We 
can then write 



{SmUp) = y^^w' pn ,{s m S„>) 



(21) 



where the weights w'i are the ideal weights given in Eqn. 
((9J and we have assumed that measurement errors are un- 
correlated with velocities. The correlation {s m s n i) can be 
calculated in the same manner as the elements of the co- 
variance matrix G (Eq. 1 19 p . 
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Once we have calculated the weights for the MV mo- 
ments, it is straightforward to calculate their correlation ma- 
trix, 



Rpq — (UpUq} — ^ ^ WpmWqn^SmSri) 



(22) 



rn n 



WpmWqnGrnn 



The moment covariance matrix can be separated into two 
parts 



(23) 



r? — »M + n( £ ) 

1 vpq — 1 l-pq ^ i ^pq j 

corresponding to the two terms in G nm . The second term 
represents the noise in the measurement of the moment 



Rpq = w pn w qn (ofj + erj) 



(24) 



The first term is due to actual motions of the objects in the 
survey and can be written as an integral over the density 
fluctuation power spectrum, 



(25) 



where the angle averaged tensor window function is 

Wp q = J w pm w qn f mn (k) (26) 

The diagonal elements W pp are the window functions for the 
individual moments u p . The window function for a partic- 
ular moment indicates which scales that moment probes. It 
can also be compared to the window function for the ideal 
moment to see the particular scales on which the moments 
differ. 

In summary, our MV method is a two-step process: for 
each moment, determine the optimal MV weights via Ea. lI5l 
then use these weights to measure the moments themselves 
fEq, H0[) . The weights depend only weakly on the power spec- 
tra and nuisance parameters (<r*), so the measured moments 
are robust. For a given moment we can also calculate its mea- 
surement variance (the diagonal elements of R p "q) and the 
cosmic variance (via Rp V q). Note that the cosmic variance 
is for the actual MV weights and not the ideal ones. The 
two variances then allow us to compare the actual moment 
amplitudes with expectations from cosmological models. 



3 DATA 

As in Paper I, the peculiar velocity data used here com- 
piles all of the major peculiar velocity su rveys published to 
date, with the exception of the survey of lLauer fc Postman! 
(|l994l ) which was found in Paper I to be inconsistent with 
other data sets. As in Paper I, we have removed outliers 
by using the predictions of the IRAS-PSCz (Point Source 
Catalogue Survey z=redshift) density f ield (for more de - 
tails see Paper I; iHudson et all |2004| ; iNeill et all |2007T >. 
Each individual survey has a characteristic MLE depth, de- 
fined as ^2r n w n /^2w n where the MLE weights are w n = 
l/(*l + * 2 «). 



This compilation, which we label "COMPOSITE", 
consists of a number of surveys, the largest of which is 
the SFlH — |- pec uliar velocity survey of spirals in the field 
and in groups dMasters et all l200rj ; ISpringob et~aH 120071 ; 
ISpringob et~aH l2009h . Here we use the data from the cor- 
rected dataset (Springob et al. 120091 ) rather than from the 



Springob et al.ll2007l ). After rejecting about 



erroneous one 

1% of the data, the sample consists of 2720 TF galaxies 
and 736 groups to make 3456 data points with characteris- 
tic depth of 35 h~ Mpc. The SFI+- 1- sample is all-sky, except 
for the Galactic plane (|fe| <> 15). As in Paper I, we find that 
all of the surveys we studied are c onsistent with each other , 
with the possible exception of the lLauer fc Postmanl l| 19941) 
(LP) BCG survey. The MV-weighted bulk flow of LP dis- 
agrees with that of the COMPOSITE catalogue on all scales. 
The level of disagreement on the larger scales, corresponds 
to 99% CL . There are independ ent reasons to question the 
LP results: IHudson et alj (I2004T ) compared, cluster by clus- 
ter, the distance to the brightest cluster galaxy derived by 
LP to that derived from the FP for other cluster members, 
and found that in a few cases, these distances differed signif- 
icantly, in the sense that the LP BCG distance was too large 
(i.e. the BCG was fainter than expected). They found that 
all of the discrepant BCGs for which HST images were avail- 
able showed strong evidence for dust. For these reasons, we 
have chosen to reject LP from the COMPOSITE catalogues. 

It is interesting to compare the results from SFI++ with 
a completely independent peculiar velocity dataset of com- 
parable depth and statistical power. The consistency be- 
tween these two catalogues, which use different data and 
methods, is an important check on our results. As in Pa- 
per I, we define the "D EEP" compilation including 103 
SNIa dTonrv et al.ll2003l ) 70 SC Tully- Fis her (TF) clusters 
(|Giovanelli et al.lll998bl : bale et al.lll999al). 56 SMAC fun- 
damental plane (FP ) clusters (IHudson et al.lll999l . l2004l ). 50 
EFAR FP clu sters l|Colless et al.ll200ll ) and 15 TF clusters 
(|Willick|[l999l ). The DEEP catalogue consists of 294 data 
points, but because these are clusters or SNe their peculiar 
velocity errors are lower per object than for SFlH — K The 
DEEP compilation covers the whole sky outside the Galac- 
tic plane and has a characteristic depth of 50 h~ 1 Mpc. 

In Paper I, we analyzed a further comp ilation, dubbed 
"SHALLOW" that consisted of the ENEAR |da Costa et"akl 
l2000l : iBernardi et~ai1l2002l ; IWegner et aT1l2003l) survey and 
a sur face brightness fluctuations (SBF) survey l|Tonrv et al.l 
l200ll ). Because for our purposes this compilation is rather 
shallow we do not analyze it separately here, but include it 
in COMPOSITE. 

In Paper I, we showed that all of these subsets (SFlH — h, 
DEEP and SHALLOW) are consistent with the same un- 
derlying velocity field, given their sparse spatial sampling. 
The formalism for comparing surveys with different window 
functions is given in Paper I (see Eq. (24) and Table 3). 

In summary, the COMPOSITE catalogue covers the 
whole sky outside the Galactic plane and has a characteris- 
tic depth of 34 h~ 1 Mpc. It is based on 4536 peculiar velocity 
measurements, making it the largest peculiar velocity cata- 
logue compiled and studied to date. 
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Figure 1. The COMPOSITE peculiar velocity catalogue on the sky (Galactic coordinates). Points are colour coded by peculiar velocity 
with red outgoing and blue infalling. While it is customary in such plots to encode the amplitude of the peculiar velocity by the size 
of the symbol, in this case symbol area is proportional to the weight. Data from SFIH — h are shown by circles, DEEP by squares and 
SHALLOW by triangles. The left panel shows the MLE weights, whereas the right panel shows the average MV weights. 



4 RESULTS 

4.1 Bulk Flow, Shear and Octupole Moments 

We calculated the MV moments for the BF, shear and oc- 
tupole components using the methods described in Sec. [2] 
above for each of the catalogues described in Sec. [3] For 
specificity, we used th e ACDM power spectrum model of 
lEisenstein fc Hul l| 19981 ) with the WMAP5 central parame- 
ters, (7 8 = 0.796, h = 0.719, fi m = 0.258 and n s = 0.963 for 
the amplitude of cosmological density fluctuations, the Hub- 
ble constant, the normalized matter density and the spectral 
index, respectively and the velocity noise it, = 150 km/s. 
The exact values of the cosmological parameters, including 
(T*, make little difference to the values of the weights. 

In principle, our MV weights allow us to match any 
choice of "ideal window", which here is assumed to have a 
Gaussian profile, parametrized by Ri. Clearly, the larger, 
denser and more geometrically complete the catalogue is 
and the smaller the velocity errors are, the closer its win- 
dow functions resemble that of the ideal window. In prac- 
tice, however, there is a compromise between two competing 
goals: one is the need to adjust the weights, as best we can, 
so that the weighted catalogue matches the ideal survey, the 
other is to keep the noise small by down-weighting objects 
with the large measurement errors. Although we have in- 
vestigated a range of scales Ri £ [10,60] /i -1 Mpc, in this 
paper, we will focus on the window with Ri — 50 /i -1 Mpc. 
This choice is a compromise between the desire to probe the 
largest possible scales, and the natural characteristic depths 
of the catalogues (~ 35 /i _1 Mpc). The observational noise 
is minimized for low Ri and becomes too large if 7?/ > 50 
/i _1 Mpc. 

Figure [1] shows the peculiar velocity data on the sky, 
both for MLE weights and for MV weights. Notice how the 
MV weights become larger in regions where the spatial sam- 
pling is poorer, such as close to the Galactic plane. The 
weighted radial distribution galaxies is shown in Figure [2] 
where it is compared to the unweighted, MLE and Ri — 50 
ft _1 Mpc ideal distributions. It is clear that the MV distri- 
bution closely approximates that of the ideal one. 

We compared the window functions (WF) (Eq. I26|) of 



all surveys and composite catalogues to the ideal survey. In 
Figures and 2] we show both the ideal and survey MV 
BF, shear and octupole squared tensor window functions 
(Eq.Hp for the COMPOSITE survey for Ri = 50/i _1 Mpc. 
The advantage of the MV moments is that they have been 
designed to be sensitive only to a narrow range of scales, 
and so we will be able to probe these scales without having 
to worry about the influence of, or aliasing from, smaller 
scales. The BF moments probes scales much larger than Ri; 
the shear responds to Ri scales and larger while the octupole 
to somewhat smaller scales (see Figures [3] and 3J) . 

When the WF's in Fig.fjare convolved by the WMAP5 
power spectrum, one finds that, for the BF statistic with 
Ri — 50/i _1 Mpc, the contribution to the integrand in Eg. 1251 
peaks at a wavenumber k ~ 0.01 h Mpc -1 , corresponding 
to wavelengths in excess of 600 /i -1 Mpc. The shear is most 
sensitive to scales where the BF window is at its minimum 
: k ~ 0.025 h Mpc" 1 , or wavelengths of 250 ft _1 Mpc. It is 
worth noting that these scales are similar to or larger than 
the very largest scales probed by Sloan Digital Sky Survey 
(SDSS) an d the 2 degree Field G alaxy Redshift Survey (2dF- 
GRS) (e.g. lPercival et al.ll2007l ). The octupole is sensitive to 
slightly smaller scales. 

Figures [3] and [4] also show that the ideal and survey 
window functions match well, especially for the lower mo- 
ments. This is also apparent in Tables [T]|3] where we give 
the values of each moment and the correlation coefficient 
{ rgrrgi ) between the ideal and their MV estimates for each 
of the moments; a correlation coefficient of unity indicates 
perfect correlations. We see that the denser the catalogue, 
the better its match with the ideal moments. 

To compare the measured moments with expectations 
from cosmology, we need to need to know what moments 
would be expected. We have only a single measurement of 
each moment (because the flow field is expanded around 
the origin, i.e. at the location of the LG) but we can calcu- 
late what values might have been measured at other loca- 
tions in a ACDM universe (assuming the Copernican princi- 
ple). For each moment we can determine its expected mean 
and variance by calculating ensemble averages over all pos- 
sible observers. Due to isotropy and homogeneity, the mean 
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Figure 2. The COMPOSITE peculiar velocity catalogue as a 
function of depth. The dashed histogram shows the distribution 
of peculiar velocity measurements. The dotted histogram shows 
the weighted histogram (renormalized to the same area) but using 
the usual MLE weights. Notice that most of the signal is driven by 
very nearby objects. The solid histogram is also weighted, but us- 
ing MV weights. The smooth curve shows the expected weighted 
radial distribution for an ideal survey i.e. oc r 2 exp I— r 2 /(2ij|)j . 
This shows that our MV weighting scheme produces the desired 
radial distribution. 



must vanish and so it is the variance that is the quantity 
of interest because it indicates the range of amplitudes for 
each moment that one would expect. We will call the corre- 
sponding standard deviation, the "cosmic root mean square" 
(CRMS) since it provides an estimate of the expected am- 
plitude of the moment for a survey with the same geometry 
and weights, but with no measurement noise. In detail, the 
CRMS is given by the diagonal elements of the covariance 
matrix R^q (eq. I25|l , which in turn depends on the weights 
(and hence the measurement noise) and on the power spec- 
trum, but does not include a contribution from the measure- 
ment noise. Tables [TJ3] include for comparison the CRMS 
of each moment around its zero mean, given the spectrum 
specified by the WMAP5 central parameters and given the 
weights as determined above. 

Turning to the values of the moments themselves, Ta- 
bles Q]|3] show that there is a remarkable agreement, within 
the errors, between all catalogues and also a good agreement 
with the CRMS (last column of the tables) for all moments 
except the Galactic y-direction component of the BF. In 
Figures and HJ] we show the bulk, shear and octupole flows 
as a function of the depth #/ for the COMPOSITE and 
SFI++ surveys. The error bars in Fig. [ 6] are for illus tration 
and are the expected rms deviation, */ ((Ui — v-i) 2 ), of the 
estimated moment from the actual value of the ideal mo- 
ment. Thus they represent a combination of measurement 



Figure 3. The window functions of the BF and shear moments 
for Rj = 50/i _:L Mpc for the COMPOSITE catalogue. The thick 
black lines are the ideal window functions for the MV compo- 
nents (since the ideal survey is isotropic, all component are the 
same) whereas the thin lines are the actual components for the 
survey. In the top panel are the BF x (short-dash), y (long-dash), 
z (dash-dot)— components. The middle panel shows the shear xx 
(short-dash), yy (long-dash), zz (dash-dot) components, the bot- 
tom panel shows the shear xy (short-dash), yz (long-dash), and 
zx (dash-dot) components. Notice how the WF virtually vanish 
at large k (small scales). All coordinates are Galactic. 



noise and deviation from the ideal window. We also show 
the expected cosmic r.m.s. (CRMS) of the flow as described 
above. As noted above, the BF is particularly sensitive to 
the large scales of the matter power spectrum, and so the 
large amplitude of the BF is suggestive of excess power on 
scales k ~ 0.01 h Mpc -1 . 

4.2 Comparison with ACDM cosmology 

To make the comparisons with cosmological models more 
precise, we compare, in a frequentist sense, the observed 
BF moments with the cosmological expectation. In other 
words, we calculate the probability that a randomly placed 
observer would have observed moments as large as the 
ones measured. Specifically, since each moment is Gaussian- 
distributed around zero mean, we calculate the x 2 for Nmom 
degrees of freedom corresponding to the number of moments, 
as given by 



%OM %OM 
2 \ ^ \ ~» 

X = 



p=l 



(27) 



where p and q specify the BF, shear and octupole compo- 
nents of the covariance matrix R pq for a specific set of values 
for the cosmological parameters we are interested in. 
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Table 1. The moment's value and its correlation coefficient ( \u\ \ u \ ^ ^ or eacn °f * ne catalogues and _R 7 = 50/i _1 Mpc ; A^mom = 19. All 
directions are given in Galactic coordinates. The last column is the expectation value of the moment given WMAP5 central parameters 
and the COMPOSITE catalogue WF's as described in the text. 





COMPOSITE 




SFI++ 




DEEP 




WMAP CRMS 


X 


86.5 ± 68.8 


0.74 


69.0 ± 95.7 


0.64 


192.7 ± 115.6 


0.51 


110.6 


y 


-404.9 ± 61.8 


0.77 


-473.6 ± 87.2 


0.67 


-320.7 ± 106.0 


0.51 


109.0 


z 


42.8 ± 37.7 


0.89 


57.7 ± 59.3 


0.80 


62.0 ± 55.8 


0.76 


105.3 


XX 


2.73 ± 1.01 


0.69 


3.36 ± 1.29 


0.62 


2.19 ± 1.76 


0.47 


1.656 


yy 


1.37 ± 0.98 


0.69 


3.72 ± 1.27 


0.63 


-0.19 ± 1.79 


0.42 


1.547 


zz 


-0.03 ± 0.68 


0.80 


2.72 ± 0.96 


0.71 


-0.72 ± 1.04 


0.67 


1.462 


x y 


0.13 ± 0.76 


0.51 


-0.71 ± 0.98 


0.42 


0.27 ± 1.29 


0.31 


0.890 


yz 


-0.95 ± 0.57 


0.63 


-1.05 ± 0.78 


0.52 


-0.71 ± 0.94 


0.40 


0.749 


zx 


1.22 ± 0.54 


0.66 


1.50 ± 0.74 


0.56 


0.98 ± 0.84 


0.47 


0.767 


XXX 


-1.2c-02 ± 2.2e-02 


0.38 


-9.3e-03 ± 2.9e-02 


0.31 


1.0e-02 ± 3.6e-02 


0.25 


1.83e-02 


yyy 


-2.4e-02 ± 1.7e-02 


0.41 


-1.9e-02 ± 2.4e-02 


0.34 


-2.2e-02 ± 2.7c-02 


0.24 


1.43e-02 


zzz 


-7.2e-03 ± l.le-02 


0.61 


-3.3e-03 ± 1.6e-02 


0.48 


-2.5e-03 ± 1.6e-02 


0.47 


1.28c-02 


xyy 


-8.2e-03 ± 1.2e-02 


0.30 


-3.3e-02 ± 1.7e-02 


0.23 


2.0e-02 ± 1.9e-02 


0.20 


8.39e-03 


yzz 


5.8e-04 ± 6.6e-03 


0.44 


-1.8e-03 ± 1.0e-02 


0.33 


8.9e-03 ± 9.6e-03 


0.30 


5.44e-03 


zxx 


7.3e-03 ± 7.8e-03 


0.45 


8.7c-03 ± l.le-02 


0.34 


-2.1e-03 ± 1.2e-02 


0.34 


6.60e-03 


xxy 


8.3e-03 ± 1.2e-02 


0.29 


5.7e-03 ± 1.6e-02 


0.24 


2.2e-02 ± 1.9e-02 


0.16 


8.24e-03 


yyz 


6.3e-04 ± 8.3e-03 


0.40 


7.7e-03 ± 1.2e-02 


0.28 


-2.5e-03 ± 1.2e-02 


0.30 


6.35e-03 


zzx 


1.2e-02 ± 7.6e-03 


0.46 


-2.5e-03 ± l.le-02 


0.35 


1.6c-02 ± l.le-02 


0.34 


6.86e-03 


xyz 


6.6e-03 ± 5.5e-03 


0.34 


9.3e-03 ± 8.2e-03 


0.25 


4.9e-03 ± 8.2e-03 


0.22 


3.72e-03 



Table 2. The same as Table Q] for N M OM = 9 





COMPOSITE 


SFI++ 




DEEP 




WMAP CRMS 


X 


101.8 ± 38.4 


0.87 


65.0 ± 46.9 


0.82 


127.6 ± 62.9 


0.70 


110.6 


y 


-362.2 ± 39.4 


0.86 


-361.6 ± 47.8 


0.81 


-326.3 ± 66.1 


0.65 


109.0 


z 


40.1 ± 30.9 


0.92 


76.6 ± 40.3 


0.87 


49.2 ± 47.8 


0.80 


105.3 


XX 


2.89 ± 0.98 


0.69 


3.94 ± 1.25 


0.63 


2.39 ± 1.61 


0.51 


1.656 


yy 


1.21 ± 0.95 


0.69 


3.90 ± 1.22 


0.64 


-0.77 ± 1.68 


0.41 


1.547 


zz 


0.10 ± 0.67 


0.80 


2.76 ± 0.93 


0.72 


-0.69 ± 1.00 


0.68 


1.462 


xy 


0.20 ± 0.74 


0.51 


-0.76 ± 0.95 


0.43 


0.22 ± 1.23 


0.31 


0.890 


yz 


-1.02 ± 0.55 


0.63 


-1.11 ± 0.76 


0.52 


-0.63 ± 0.86 


0.41 


0.749 


zx 


1.44 ± 0.53 


0.67 


1.48 ± 0.72 


0.56 


1.32 ± 0.80 


0.48 


0.767 



The results are presented in Table [4] where we show, 
for each of the catalogues, the percent probability of get- 
ting a larger x 2 , P{> X 2 )- I n the table, we break down the 
probabilities for Nmom=3 (BF alone), Nmom=9 (total, BF 
and shear) and for Nmom=19 into total, BF, shear and oc- 



tupole. This has the benefit of showing us clearly which part 
of the flow is in agreement or disagreement with standard 
cosmology. Table [3] shows that although the BF is highly 
unlikely, both the shear and octupole moments amplitudes 
are more or less what is expected. For Nmom=19, the total 
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Table 3. The same as Table [T] for A^mom = 3 





COMPOSITE 


SFI++ 




DEEP 




WMAP CRMS 


X 


97.4 ± 38.3 


0.87 


71.7 ± 46.4 


0.82 


120.0 ± 61.2 


0.71 


110.6 


y 


-361.9 ± 39.2 


0.86 


-363.3 ± 46.9 


0.82 


-333.8 ± 65.4 


0.65 


109.0 


z 


41.2 ± 30.7 


0.92 


73.1 ± 39.8 


0.87 


44.3 ± 47.1 


0.80 


105.3 




Figure 4. The same as \3\ for the octupole moments. The thick 
black lines represent the ideal WF. In the top left panel are the 
xxx (short-dash), yyy (long-dash), zzz (dash-dot); the top right 
panel shows the xxy (short-dash), yyz (long-dash), zzx (dash- 
dot), the bottom left panel shows the xyy (short-dash), yzz (long- 
dash), and zxx (dash-dot); the bottom right panel shows the xyz 
(short-dash) component of the WF's. 

probability of getting the observed \ 2 or higher is ~ 20%, 
the probability of the shear alone is € [10, 50]% and the 
octupole ~ 80%, whereas getting the BF is <> 0.5%. 

Note that the BF is discrepant at the 99.78% and 
97% levels for the two independent subsamples (SFI++ and 
DEEP respectively). Recall that the SFI++ peculiar ve- 
locities are based on the Tully-Fisher relation applied to 
field and group spirals, whereas the DEEP compilation is 
based on a variety of other methods, particularly Type la su- 
pernovae and the Fundamental Plane relation in early-type 
galaxies. The agreement between these independent datasets 
indicates that the BF is not a result of systematic errors. 

Thus when considered by itself, the BF disagrees with 
ACDM at > 98% CL, but there is no disagreement when all 
19 moments are considered together. This is a consequence 
of the fact that we have only one component, the Galactic 
y-direction of the BF, that is much higher than expected. 
When considering only three components of the BF, then 
clearly if one disagrees, the \ 2 wm be much higher (probabil- 



Figure 5. The BF and shear moments of the COMPOSITE (red 
solid line) and SFIH — (- (blue dashed line) catalogues as a function 
of Rj, all in Galactic coordinates. Note that the data points are 
not independent. The top panels are the BF velocities to the 
Galactic x (left), y (center) and z (left) directions. The error bars 
are as described in the text. In the middle and bottom panels we 
show the shear components in the Galactic directions xx (middle 
left), yy (middle center) and zz (middle right) and the xy (bottom 
left) yz (bottom center and zx (bottom right). The solid black 
line, without the error bars, in each panel is the WMAP5 ACDM 
cosmic rms (CRMS). 



ity much lower) than if only one out of nineteen amplitudes 
is in conflict. This trend is clearly shown in Table [4] 

On the other hand, as noted above, for the shear and oc- 
tupole moments, the measurement uncertainties are greater 
than for the BF because these moments are essentially 
derivatives of the (noisy) flow field. Thus the shear, for ex- 
ample, is not as powerful a cosmological probe as the BF, at 
least with current datasets. One way to quantify this is to 
compare the WMAP5 CRMS to the measurement noise (see 
Figs. EJ M and Tables [T] [2} . In the case of the BF, this ratio 
is in the range 1.5 — 2.5, whereas for the shear it's 1 — 1.2 and 
for the octupole the values are only 0.5 — 1. Thus while we 
have shown that these higher moments are not much differ- 
ent than expected, due to the noise we do not have much of 
a handle on their actual values. Consequently, the discrim- 
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This BF value also agrees remark ably wel l with 
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Figure 6. The same as Fig.[5]for the octupole moments. Galactic 
xxx (top left), yyy (top center) and zzz (top left), xxy (middle 
left), yyz (middle center) and zzx (middle right), the xyy (bottom 
left) yzz (bottom center and xyz (bottom right). 



inatory power of these higher order moments (and their x 2 
values) is less than for the BF, since the contribution to the 
X 2 from measurement-noise-dominated moments is always 
going to be about unity. 

Finally, note that the shear tensor tabulated in Tables [T] 
and [2] and shown in Figure [5] has 6 components. In practice, 
however, the trace of the shear (i.e. the average expansion 
rate) is typically not measurable by peculiar velocity sur- 
veys, but is a free parameter that is usually adjusted to be 
close to zero by requiring that there be no net inflow or out- 
flow of peculiar velocity tracers. Therefore, it is more conser- 
vative to consider the traceless components of the shear. If 
we subtract off 1/3 of the trace from each component on the 
diagonal (xx, yy and zz), we see that the measured shears 
are similar to the WMAP5 expected rms values. 



5 DISCUSSION 

Our total BF, calculated for the COMPOSITE catalogue 
at Ri = 50 ft^Mpc, with N MO m=19 gives \v\ = 416 ± 78 
km/s towards Galactic I = 282° ± 11° and b = 6° ±6° which 
is in disagreement w ith the expectations of the WMAP5 
jPunklev et al.ll2009l ) cosmology at the 98-99.5% CL. This 
result, however, is in excellent agreement with the results 
found previously (Paper I) estimating only the BF. The 
small magnitude of the shear and octupole moments sug- 
gests that our previous result was not due to, for example, 
aliased small-scale power contaminating the BF measure- 
ment. This shows that our orthogonality procedure (Eq. [8]) 
works very well. 



other peculiar velocity estimates (e.g. Hudson 19941; 



2005 



Sarkar et al 



Hoffman et"al] 1200 1| ; IZaroubi et al.l l200ll; Pike fc Hudson! 



20071 ; iFeldman fc Watkins! I200I 



Lavaux et al l l201Ch . though these estimates were derived 



with different catalogues, methodologies and assumptions. 

The results for the shear and octupole are consistent 
with the hypothesis that the power is not unusual on scales 
smaller than the very large ones probed by the BF. As men- 
tioned above, another way to mak e such a compar i son is 
via supercluster infall. For example, iPike fc Hudson! (|2005h 
found (fi m /0.3)°- 55 (78 = 0.80 ± 0.05, consistent with the 
mean WMAP5+BAO+SN values (Q m /0.3)°- 55 o - 8 = .77 ± 
0.035 jKomatsu et al.ll2009l ). [Abate fc Erdogdul |2009l ). who 
studied the velocity correlation function of SFI++, a statis- 
tic that mixes a range of scales, also found consistency with 
WMAP5 cosmological parameters. Thus it does seem to be 
the BF which is unusual. 

IKashlinskv et all (|2008l . 



Recently. 



lAtrio-Barandela et all (|20ld ) 



l201Ch and 



claimed to have detected 
a dipole in filtered WMAP CMB temperature maps mea- 
sured at the locations of rich clusters. The magnitude of 
this cluster-temperature dipole is 2.8 ± 0.7 /iK. The authors 
interpret this as being due to the kinetic Sunyaev-Zel'dovich 
(kSZ) effect, implying a large scale flow with a bulk velocity 
in the range € [600, 1000] km/s towards I = 283° ± 14°, 
b = 12° ± 14°. Although the scale of their cluster catalogues 
(z ~ 0.1 — 0.2) is of order three or more times the scales 
that our catalogues probe, this kSZ flow is in excellent 
directional agreement with the one we detected here and in 
Paper I. The magnitude is larger, but the authors caution 
that the magnitude is systematically uncertain. We note 
that the BF that we detect does not seem to level off yet 
and there is no clear sign of convergence to the CMB frame 
in our catalogues (see Figure [5J top center panel). We by 
no means claim to verify the IKashlinskv et al. | (|2008l . l2010l ) 
results, however, we certainly can not refute it with our 
data. 

The motion we detect is not due to nearby sources, such 
as the Great Attractor (distance of ~ 40ft. Mpc), but rather 
to sources at greater depths that have yet to be fully iden- 
tified. The largest known mass concentration, the Shapley 
supercluster, does not see m to be massive eno ugh to cause 
a flow of this magnitude (|Ravchaudhurvlll989l ). It is likely 
that the flow arises both from various mass concentrations 
in the Galactic y-direction as well as under-dense regions in 
the opposite direction. Currently, there is no peculiar veloc- 
ity survey in existence that is deep enough to resolve the 
source(s) of the flow, if indeed it is a potential flow. 

A good measure for the distance scale to the sources re- 
sponsible for mos t of the BF is the ratio of the BF to shear 
(|Lilie et al.l Il98rj ; iKaiserl fl99ll ). This characteristic depth 
Di oc Ui/ua with constant of proportionality of order one 
should be valid even if multiple sources, including under- 
dense regions, as expected for a Gaussian random field, are 
responsible for the flow. In Fig. [7] we show the effective dis- 
tance D y oc Uy/uyy to an attractor for the COMPOSITE 
survey as a function of Ri. Nearby, where Ri <J 30 ft _1 Mpc, 
the moments are dominated by the nearby velocity field and 
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Table 4. The total observed P(> X 2 ) m percent for Nmom= 3, 9 and 19 and for the BF, shear and octupole moments for each, for 
Rj = 50/i -1 Mpc, and the WMAP5 central parameters U m = 0.258 <r 8 = 0.796. 





Wmom = 3 
BF 


Total 


JVmom = 
BF 


9 

shear 


Total 


-?Vmo 
BF 


M — 19 
shear 


octupole 


COMPOSITE 


1.89 


6.01 


1.81 


41.76 


17.00 


0.50 


52.60 


78.33 


SFI++ 


3.11 


1.73 


3.22 


7.70 


16.19 


0.22 


11.22 


89.38 


DEEP 


6.02 


30.41 


6.29 


82.62 


55.54 


3.18 


91.22 


81.61 
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400 - 
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Figure 7. The expected distance to the dominant gravita- 
tional sources (solid line) (u y /u yy ) as a function of Rj for the 
COMPOSITE survey as well as an expected distance using the 
WMAP5 central parameters (solid line, no errorbars) discussed 
in the paper. For comparison we also show u x /u xx (dashed line) 
that agrees with the WMAP5 expectations within errors. 



show the effect of the Great Attractor at ~ 35 /i _1 Mpc. 
Then there is a sudden rise around 30/i -1 Mpc that sug- 
gests that on scales beyond the Great Attractor, the flow 
is driven by even larger and more distant mass concentra- 
tion^) in the same general direction. The uncertainties are 
large, however, at Ri = 50/i~ 1 Mpc we find a characteristic 
distance of 296 ±62 h~ 1 Mpc, much larger than the expected 
WMAP5 value of 82 h^Mpc. For comparison we show the 
estimated u x /u xx which is consistent with the WMAP5 ex- 
pected results. 

Note, however, as discussed above, it is sensible to sub- 
tract 1/3 of the trace of the shear tensor from each diagonal 
component. In this case, the yy component of the shear ten- 
sor is consistent with zero, implying that the attractor is at 
infinity. Thus we cannot exclude the cohe r ent flow on much 
larger scales claimed bv lKashlinskv et al l (|2008l . l20TOl ). 



6 CONCLUSIONS 

We have estimated the MV BF, shear and octupole moments 
and showed that they are minimally sensitive to aliasing 
from small-scale power. We have applied the MV formal- 
ism to a number of compilations of recent peculiar velocity 
surveys. The MV window functions are similar to the ideal 
window functions and our octupole moments orthogonaliza- 
tion procedure works well in a sense that including them 
does not affect the BF much. 

The various peculiar velocity surveys are consistent 
with each other and all show a large BF which does not 
follow expectations from the WMAP5-normalized ACDM 
model. Specifically, we have shown that the BF within a 
Gaussian window of radius Ri = 50fo _1 Mpc has a magni- 
tude of 416 ± 78 km/s towards Galactic I = 282° ± 11° and 
b = 6 ° ± 6° in disagreement with WMAP5 (jPunklev et all 
I2009T ) at the 99.5% CL. This flow is consistent with be- 
ing relatively "cold" (both shear and octupole moments are 
in agreement with expectations) and is not due to nearby 
sources. If we include all nineteen octupole (nine shear) mo- 
ments, we disagree with the WMAP5 expectations at only 
the 83% (94%) CL, respectively. We also found that the 
sources responsible for the BF are at an effective distance 
of > 300 ft, - Mpc, too far to identify in existing all-sky red- 
shifts surveys. 

If the flow is a potential flow, there are various possible 
explanations as to its source. While it is always possible that 
the result is due to systematic errors in the data, we tend to 
discount this possibility since the BF is seen in many inde- 
pendent surveys with various distance indicators, method- 
ologies and geometries. It may be that we live in a ACDM 
Universe but that we happen to live in one of the very rare 
volumes that exhibits this flow (<J 1%), or we may live 
in a Universe with more large-scale power than WMAP5- 
normalized ACDM. In either of these last two cases, we may 
be able to find the source or sources of the flow in future 
redshift surveys. It is of course likely that the sources are a 
combination of over- and under-dense regions and not just 
a single large mass concentration. 
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